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Abstract 

Background: Genetic interactions within hybrids influence their overall fitness. Understanding the details of these 
interactions can improve our understanding of speciation. One experimental approach is to investigate deviations 
from Mendelian expectations (segregation distortion) in the inheritance of mapped genetic markers. In this study, 
we used the copepod Tigriopus californicus, a species which exhibits high genetic divergence between populations 
and a general pattern of reduced fitness in F2 interpopulation hybrids. Previous studies have implicated both 
nuclear-cytoplasmic and nuclear-nuclear interactions in causing this fitness reduction. We identified and mapped 
population-diagnostic single nucleotide polymorphisms (SNPs) and used these to examine segregation distortion 
across the genome within F2 hybrids. 

Results: We generated a linkage map which included 45 newly elucidated SNPs and 8 population-diagnostic 
microsatellites used in previous studies. The map, the first available for the Copepoda, was estimated to cover 75% 
of the genome and included markers on all 12 1 californicus chromosomes. We observed little segregation 
distortion in newly hatched F2 hybrid larvae (fewer than 10% of markers at p < 0.05), but strikingly higher 
distortion in F2 hybrid adult males (45% of markers at p < 0.05). Hence, segregation distortion was primarily caused 
by selection against particular genetic combinations which acted between hatching and maturity. Distorted 
markers were not distributed randomly across the genome but clustered on particular chromosomes. In contrast to 
other studies in this species we found little evidence for cytonuclear coadaptation. Instead, different linkage groups 
exhibited markedly different patterns of distortion, which appear to have been influenced by nuclear-nuclear 
epistatic interactions and may also reflect genetic load carried within the parental lines. 

Conclusion: Adult male F2 hybrids between two populations of T. californius exhibit dramatic segregation 
distortion across the genome. Distorted loci are clustered within specific linkage groups, and the direction of 
distortion differs between chromosomes. This segregation distortion is due to selection acting between hatching 
and adulthood. 



Background 

One way in which the integrity of species can be main- 
tained is by intrinsic postzygotic barriers which reduce 
the fitness of hybrid offspring. Understanding how 
genetic interactions within interpopulation hybrids influ- 
ence fitness is hence an important step towards under- 
standing how new species arise. When gene flow 
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between populations is restricted, they can diverge as a 
result of selection and drift. One outcome of this diver- 
gence may be the accumulation of 'Dobzhansky-Muller 
incompatibilities', mutations which only become deleter- 
ious when placed on a novel genetic background in 
hybrids between the populations [1,2]. As first genera- 
tion (Fl) hybrids contain a full haploid nuclear genome 
from each of the parental populations, the negative fit- 
ness consequences of Dobzhansky-Muller incompatibil- 
ities may not be expressed until the second hybrid (F2 
and backcross) generations, when co-adapted blocks of 
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the genome have been broken up by recombination. 
Studies have now identified regions of the genome, and 
in some cases specific genes, involved in such deleter- 
ious Dobzhansky-Muller interactions in hybrids [3-6]. In 
contrast, some authors have observed evidence for 
favorable interactions between divergent parental gen- 
omes in interpopulation and inter-species hybrids 
[7-10], suggesting that recombination of parental gen- 
omes also has the potential to increase fitness. Examina- 
tion of the patterns of segregation distortion at genetic 
markers distributed across the genome within hybrids 
can both help identify candidate regions for Dobz- 
hansky-Muller incompatibilities [11], and provide 
insights into other types of genetic influence on hybrid 
fitness. 

The harpacticoid copepod Tigriopus californicus, an 
easily cultivated species that is common in supra-littoral 
splash pools on rocky outcrops along the west coast of 
North America, has become an established model for 
investigating the genetic basis of hybrid breakdown. 
There is little gene flow between populations, even 
those separated by relatively short distances [12,13]. 
Divergence in mtDNA nucleotide sequence between 
populations separated by less than 60 km can exceed 
20%, which encompasses both synonymous and nonsy- 
nonymous sites and reflects both this low gene flow and 
an unusually high mtDNA substitution rate compared 
to nuclear genes [14,15]. The level of divergence at 
nuclear loci is still substantial, but is between 5.8 and 
38.2-fold lower than mitochondrial divergence between 
the same geographic locations [14]. Despite such 
extreme genetic differentiation, almost all populations 
can easily be crossed to produce viable hybrids [16,17]. 
Interpopulation crosses typically exhibit a pattern of 
slightly increased fitness in the Fl generation, followed 
by lowered fitness in the F2 and later hybrid genera- 
tions, as demonstrated by a variety of measurements 
including fecundity, development time and survivorship 
[18-20], response to osmotic stress [21,22], cytochrome 
oxidase activity [23] and mitochondrial ATP production 
[22]. This pattern suggests a role for Dobzhansky-Muller 
incompatibilities in causing the hybrid breakdown. The 
degree of fitness reduction in F2 hybrids shows a signifi- 
cant association with genetic and geographic distance 
between populations [16]. 

Much work has focused on the role of nuclear-mito- 
chondrial incompatibilities in this hybrid breakdown. 
Ellison and Burton [20] found clear evidence of this 
intergenomic interaction. They backcrossed F3 hybrid 
females, which exhibited lowered fitness, to pure males 
from each parental population: they found that restoring 
a full nuclear haplotype in this way restored fitness only 
if the mtDNA from the matching population was also 
present. At least two sets of candidate loci for nuclear- 



mitochondrial interactions have attracted investigation. 
First, there has been much interest in the interaction 
between nuclear and mitochondrial encoded proteins 
involved in the oxidative phosphorylation (OXPHOS) 
pathway. In [24], the authors observed a general trend 
towards loss of fitness, associated with reduced ATP 
production, in hybrid lines compared to parental con- 
trols. They found that those OXPHOS enzyme com- 
plexes that contained both mitochondrial and nuclear 
components exhibited reduced activity in hybrid lines, 
whilst those containing only nuclear components did 
not, implicating nuclear-mitochondrial interactions in 
generating the observed results. Several studies have 
found strong evidence for nuclear-mitochondrial coa- 
daptation affecting the activity of cytochrome c oxidase 
(COX), an enzyme with both nuclear and mitochondrial 
encoded units, and its interaction with the nuclear 
encoded cytochrome c [23,25,26]. Population specific 
interactions between COX and cytochrome c have been 
traced to the level of individual amino acid substitutions 
[27]. Second, Ellison and Burton [20] have suggested 
that there is also coadaptation in the mitochondrial 
transcription apparatus. Mitochondrial RNA polymerase 
(mtRPOL, nuclear-encoded) initiates mtDNA transcrip- 
tion by binding specific mtDNA sequence in the non- 
coding "control region." Using inbred hybrid lines, the 
authors [22] found that lines with mtRPOL and mtDNA 
from the same population exhibited parental OXPHOS 
expression patterns, whilst those with mismatching 
mtRPOL and mtDNA did not. 

However, cytonuclear coadaptation is clearly not com- 
plete within all T. californicus populations [23,28] and 
evidence suggests that interactions within the nuclear 
genome also play a role in hybrid breakdown. For exam- 
ple, Ellison and Burton [22] found that, for hybrid lines 
containing mtDNA from a San Diego population (SD), 
parental expression was not restored with the matching 
mtRPOL, suggesting that at least one additional, epista- 
tically interacting nuclear locus is present in this popu- 
lation but not in the experimental lines. Willett and 
Berkowitz [29] found dramatic segregation distortion for 
two nuclear encoded homologues of malic enzyme in F2 
hybrids, which could be attributed to selection against 
particular genotypes acting between hatching and sexual 
maturity. Willett [30] examining segregation patterns of 
three nuclear-encoded OXPHOS genes, found evidence 
for complex nuclear epistatic interactions, but not 
nuclear-mitochondrial interactions, affecting survival to 
adulthood of F2 hybrids. Edmands and colleagues 
[28,31] similarly implicated epistatic interactions 
between different chromosomes in influencing segrega- 
tion distortion in backcross and F2 hybrids. 

In this study, we investigated the genetic composition 
of T. californicus interpopulation hybrids in more detail. 
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We made crosses between two populations, SD and SC, 
that are 640 km apart and more than 20% divergent in 
mitochondrial COI and CYTB sequence [15,16]. These 
populations display the classic pattern of decreased sur- 
vivorship of F2 interpopulation hybrids compared to 
parentals (Edmands and coworkers, unpublished data), 
and have been used in previous studies. We examined 
segregation pattern of eight microsatellite markers and 
45 newly elucidated single nucleotide polymorphism 
(SNP) markers, distributed over the genome, in the F2 
individuals. We compared patterns of segregation distor- 
tion between newly hatched nauplii (larvae) and adults 
in order to distinguish distortion caused by differential 
post-hatching viability from distortion caused by either 
meiotic drive, differential gametic fertilization success or 
pre-hatching mortality. 

Methods 

SNP discovery and genotyping 

A cDNA library was constructed from RNA extracted 
from a mass copepod sample collected from the San 
Diego population (SD, 32° 45'N, 117° 15'W) of T. cali- 
fornicus. This sample (approximately 1 g wet weight) 
included adults of both sexes as well as eggs, nauplii 
and copepodid developmental stages. Total RNA was 
extracted using Tri Reagent (Sigma Chemical) using the 
manufacturer's protocol. Purification of mRNA from 
total RNA was achieved using a Qiagen mRNA isolation 
kit. Approximately 6 ug of mRNA was then used to cre- 
ate the cDNA library using the Zap Express cDNA 
Synthesis Kit (Stratagene), following manufacturer's pro- 
tocols. Cloned cDNAs were recovered by in vivo exci- 
sion of the pBluescript* phagemid in E. coli, and 1100 
positive clones (based on blue/white screening) were 
randomly selected for sequencing. Inserts from the posi- 
tive clones were PCR amplified using M13 primers and 
then sequenced using T3 and T7 primers. Sequencing 
was carried out on an Amersham MegaBACE 500 
sequencer with Amersham's ET Dye-Terminator 
chemistry. 

From the cDNA sequences, we selected non-mito- 
chondrial sequences over 700 bp in length and designed 
primers to amplify 500-700 bp portions of these 
sequences using Primer 3 (http://primer3.sourceforge. 
net/). For all primer design, we specified an optimal pri- 
mer length of 20 bp and optimal melting temperature of 
60°C. We also designed primers to amplify three por- 
tions of the mitochondrial genome, using mitochondrial 
sequences for T. californicus deposited in GenBank. Pri- 
mers were manufactured by Operon Biotechnologies 
(http://www.operon.com) and IDTDNA (http://www. 
idtdna.com). We tested each primer pair on two indivi- 
duals from each of three populations: SD, SC (Santa 
Cruz, California, 36°57'N, 122°03'W), and PBJ (Punta 



Baja, Baja California, 36°57'N, 122°03'W, used in a sepa- 
rate study). DNA was extracted by placing copepods in 
50 ul lysis buffer (10 mM Tris pH 8.3, 50 mM KC1, 
0.5% Tween 20) with 200 ug/ml Proteinase K and incu- 
bating at 65°C for 2 hr followed by 100°C for 15 min. 
We amplified products in 25 ul reactions using the fol- 
lowing recipe: 2.5 ul template DNA; 2.5 ul each forward 
and reverse primers (10 mM each); 2.5 ul premixed 
dNTPs (2.0 mM each); 2.0 ul MgCl 2 (25 mM); 2.5 ul 10 
X buffer, 10.4 ul DNA-free water, 0.1 ul Taq polymer- 
ase, and the following reaction conditions: a denatura- 
tion step of 94°C for 5 minutes; 35 cycles of 94°C for 30 
sec, 55°C for 30 sec, 72°C for 30 sec (45 sec for mtDNA 
primers), followed by a final extension step of 72°C for 
5 minutes. PCR products were visualized on 1.8% agar- 
ose gels with ethidium bromide staining. Where a pri- 
mer pair successfully produced a single amplification 
band in SD and at least one of the other test popula- 
tions, the product was sequenced. Products were sent to 
the High-Throughput Genomics Unit at the University 
of Washington (http://www.htseq.org) for exo/sap clean- 
up followed by unidirectional sequencing. Two indivi- 
duals per population were considered sufficient for pre- 
liminary identification of population-specific SNPs as 
geographically isolated T. californicus populations tend 
to be genetically homogeneous [32]. Levels of poly- 
morphism were subsequently determined by genotyping 
at least 28 individuals from each of the three 
populations. 

We used Sequencher v. 4.6 (http://www.genecodes. 
com) to align sequences and searched for SNPs diagnos- 
tic between both SD and SC, and SD and PBJ, with >25 
bp of invariant sequence either upstream or downstream 
to enable the design of single base extension primers. In 
addition, we screened T. californicus mtRPOL sequence, 
obtained from GenBank, to identify suitable SNPs 
between SD and SC and SD and PBJ. We selected one 
suitable SNP site for each separate nuclear sequence 
and three from the mitochondrial sequences. SNP and 
flanking sequence information were provided to Jeffrey 
Conroy, Genomic Shared Resources, Roswell Park Can- 
cer Institute, Buffalo, NY, who performed primer design 
using iPLEX Gold software (Sequenom, San Diego, CA). 

All SNP genotyping was performed at Roswell Park 
Cancer Institute, using the iPLEX Gold Assay on a Mas- 
sARRAY Compact (Sequenom). We provided unpuri- 
fied, dried-down, lysis extracts (10 ul for adults, 20 ul 
for nauplii). Diagnostic utility of SNPs was confirmed by 
genotyping a larger sample of copepods from SD, SC 
and PBJ. Accuracy of SNP calls in heterozygotes was 
tested by including a blind set of known Fl SD x SC 
and SD x PBJ individuals. We also investigated the uti- 
lity of the SNPs for T. californicus research in general 
by genotyping them for 13 additional populations, 
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ranging from Washington to Baja California (Additional 
Files), the majority of which have been used in previous 
studies [17,32,33] 

Experimental Crosses 

Reproductive biology of T. californicus is well estab- 
lished [33-35]. Males guard immature females by clasp- 
ing them with their antennae until the female completes 
her terminal molt; the female is then inseminated and 
released. Unmated females are therefore easily obtained 
by separating clasped pairs. Females mate only once and 
use stored sperm to fertilize multiple clutches of eggs, 
with each female producing an average of ~300 progeny 
[35]. Inbred 'isofemale' lines are therefore established by 
isolating a single fertilized female and allowing offspring, 
including those from overlapping generations, to mate 
freely. 

Mapping crosses were initiated using isofemale lines 
from populations SC and SD (one line per population) 
that had been inbred for 9-10 months (minimum gen- 
eration time is 23 days at 20°C [36]). All crosses were 
performed in 60 x 15 mm Petri dishes containing 10 ml 
growth medium (1 liter seawater filtered through a 37 
um filter, 0.1 g ground Tetramin fish food, 0.1 g pow- 
dered Spirulina). One virgin female from the SD line 
and one adult male from the SC line were placed in 
each Petri dish. The male was removed when the first 
eggsac was observed and the female moved to a new 
Petri dish after each eggsac hatched. The Fl offspring 
from each clutch were allowed to grow to maturity in 
the dish. When clasped pairs were observed, they were 
individually moved to a fresh Petri dish, where they pro- 
duced the first batch of the F2 generation. Males were 
removed once an eggsac was observed and females 
again moved to a new dish after each eggsac hatched. 
All dishes were housed in a 20°C incubator with a 12 h 
light: 12 h dark cycle. Extra Spirulina was added when 
supplemental food was considered necessary. As crosses 
were initiated using SD females, all F2 and backcross 
individuals contained the SD mitochondrial haplotype. 

To distinguish between segregation distortion caused 
by meiotic drive or differential gametic fertilization suc- 
cess, and that caused by post-zygotic selection against 
particular genetic combinations, we compared results 
from newly-hatched nauplii to those from adult males. 
Newly hatched F2 nauplii for genetic analysis were 
obtained by removing a late-stage (orange-colored) egg- 
sac from a female and allowing nauplii to emerge in a 
drop of seawater in a Petri dish. Nauplii were then killed 
by flooding the dish with lysis buffer and individual nau- 
plii were transferred in 2-5 ul of liquid to a 200 ul PCR 
tube containing 20 ul of lysis buffer with 200 ug/ml Pro- 
teinase K. Tubes were incubated at 65°C for 1 hr fol- 
lowed by 100°C for 15 min and extracts were then 



stored by freezing at -70°C. Adult males were obtained 
by allowing F2 families to reach maturity in the Petri 
dish. All males in a dish were rinsed in deionized water 
and frozen whole in 200 ul PCR tubes. DNA was subse- 
quently extracted in 50 ul of lysis buffer with 200 ug/ml 
Proteinase K using the protocol described above. Geno- 
typed nauplii (n = 190) were the offspring of four differ- 
ent Fl pairs; genotyped adult males (n = 205) were the 
offspring of 25 different Fl pairs. 

Resources limited genotyping to only one of the two 
reciprocal crosses and to adults of only one gender 
(males). Males were chosen to eliminate the possibility 
of amplifying sperm or zygote alleles in fertilized 
females. Tigriopus californicus lacks heteromorphic sex 
chromosomes [37,38] and sex determination, although 
still poorly understood, may involve both additive 
genetic and environmental components [39,40]. Previous 
studies of this species have found differences between 
the sexes in segregation distortion [25,30,31,41], most 
commonly finding less distortion in males; however 
these differences are not systematic, varying with the 
populations used, direction of cross, locus and replicate. 
We acknowledge that focusing only on adult males lim- 
its the scope of data interpretation; as the first genome- 
wide assessment of larval and adult F2 genotype fre- 
quencies in this system, our study nevertheless offers 
important insight into the mechanism driving segrega- 
tion distortion. 

As female T. californicus do not undergo chromoso- 
mal recombination [37], we were also able to produce 
non-recombinant backcrosses (NR-BC) to confirm that 
identified linkage groups were on different chromo- 
somes. Fl females were backcrossed to SC males. As 
before, NR-BC families were allowed to develop to 
maturity in Petri dishes and adult males (n = 39, the off- 
spring of seven different pairs) frozen prior to genetic 
analysis. 

Microsatellite genotyping 

In addition to the SNPs, we genotyped F2 and NR-BC 
individuals for nine microsatellite loci (30, 197, 228, 480, 
558, 1202, 1203, 1555, 56J2 [42]) that exhibited fixed 
size differences between the SD and SC populations in a 
test sample (n = 24 each). These markers enabled SNP 
linkage groups to be anchored to eight of the 12 T. cali- 
fornicus chromosomes numbered by Harrison and 
Edmands [31]. Microsatellites were amplified individu- 
ally in 12 ul reactions containing: 0.5 ul template DNA; 
0.3 ul fluorescently labeled forward primer (10 mM); 1 
ul reverse primer (10 mM); 1.2 ul premixed dNTPs (2.0 
mM each); 1.2 ul MgCl 2 (25 mM); 1.2 ul 10X buffer, 
6.04 ul DNA-free water, 0.06 ul Taq polymerase, with 
the Thermocycler conditions as described for SNP dis- 
covery above (annealing temperature for 480 = 61°C, for 
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all other microsatellite loci 55°C). PCR products were 
multiplexed into two pools and run on a CEQ 8000 
capillary sequencer (Beckman Coulter, Fullerton, CA, 
USA). Allele sizes were scored manually. 

Statistical Analysis 

Separate linkage maps for adults and nauplii were gen- 
erated using Map Manager QTX, using the Kosambi 
map function with linkage criterion set at p < 0.001. We 
estimated corrected map length (L) and coverage (c) 
using the methods of [43] and [44] as follows. For each 
linkage group, we added two cM to the total length and 
then multiplied this by (m+l)/(m-l), where m is the 
number of markers on each group. We then estimated 
genomic coverage as c = 1 - e 2 dn/L where d = mean 
intermarker distance and n = total number of markers 
assigned to linkage groups. To examine segregation dis- 
tortion in nauplii and adults, we compared observed sin- 
gle-locus genotype frequencies to those expected from 
Hardy- Weinberg equilibrium using the x -squared good- 
ness-of-fit test. Additionally, we investigated whether 
pairwise epistatic interactions were contributing to seg- 
regation distortion by looking at the co-occurrence of 
genotypic classes at different SNP loci. For each pair of 
physically unlinked loci we calculated expected frequen- 
cies of each genotypic combination (SDSD/SDSD; 
SDSD/SDSC; SDSD/SCSC; SDSC/SDSC; SDSC/SCSC; 
SCSC/SCSC) from the observed frequency of each geno- 
type at each locus. We then compared observed and 
expected frequencies of genotypic combinations using a 
X-squared goodness-of-fit test. We also used ^-squared 
tests to examine chromosome-level segregation distor- 
tion and pairwise linkage disequilibrium between chro- 
mosomes in the small non-recombinant backcross 
sample. 

Results 

We generated over 20,000 bp of sequence data for each 
of the three T. californicus populations included in the 
SNP discovery step. Mean observed nuclear sequence 
divergence between SD and SC, and between SD and 
PBJ, was 3.8%, and 3.5% respectively. Sequence data was 
submitted to NCBI GenBank (Additional file 1, Table 
SI). We designed assays for 51 nuclear SNPs and three 
mitochondrial SNPs (Additional file 1, Table SI), which 
could be multiplexed into two pools for iPLEX Gold 
genotyping. Of these, 49 SNPs were identified as diag- 
nostic between SD and SC. On the basis of results from 
Fl individuals, we rejected a further four SNPs from 
this study because SD-SC heterozygotes were under- 
called, leaving 45 nuclear markers (Additional file 1, 
Table 1). The majority of SNPs could be scored in all 
but one of 14 additional T. californicus populations, and 
generally appeared fixed in each population (Additional 



file 1, Table SI). However, we note that most of our 
assays were optimized for just two alternative SNPs, 
meaning that the presence of additional undetected SNP 
alleles that did not occur in our three original popula- 
tions cannot be ruled out: 2% of the total non-indel 
polymorphic sites observed in our sequence data exhib- 
ited three alleles. Half of all SNPs failed in PA, the only 
population known to exhibit reproductive isolation from 
others [17]. 

One microsatellite (1202) was discovered to be non- 
diagnostic between SD and SC and was not included in 
the linkage map. Parents of the non-recombinant back- 
cross individuals carried distinct alleles at this locus and 
therefore we were still able to use it to anchor one link- 
age group to a named chromosome. 

For both nauplii and adults, nuclear SNPs and micro- 
satellites formed 11 linkage groups and one unlinked 
locus, corresponding to the 12 chromosomes of T. cali- 
fornicus (Figure 1) [37]. All linkage relationships were 
strongly supported, with all LOD scores greater than 11 
and most much higher (Additional file 2, Table S2). 
However LOD scores between microsatellites markers 
and SNP markers were generally lower than between 
SNP pairs, reflecting the fact that fewer individuals were 
successfully genotyped for microsatellite loci (mean total 
n for SNP loci = 389; for microsatellite loci = 298). 
Results from non-recombinant backcross individuals 
confirmed each linkage group to be on a separate chro- 
mosome. Inclusion of the microsatellites enabled eight 
chromosomes named by Harrison and Edmands [31] to 
be identified (Figure 1). Total length of the nauplii link- 
age map, calculated by summing all inter-marker dis- 
tances, was 266.7 cM Kosambi, with a mean between- 
marker distance of 6.5 cM and a maximum distance of 
21.2 cM. Total adult map distance was very similar at 
263.4 cM, however length differences were present 
between homologous linkage groups in the two maps, 
reflecting increased segregation distortion in the adults. 
Most markers were co-linear between the two maps. 
The greatest variation was caused by the shifting of two 
microsatellite markers, which may be an artifact of the 
smaller sample size for this marker type. Maps gener- 
ated using the SNPs alone (Additional file 2, Table S2) 
were highly congruent, with the exception that marker 
TC085 was not assigned to a linkage group in the nau- 
plii. Including the microsatellites, corrected map length 
for the nauplii was 484.8 cM, and estimated genomic 
coverage for this map was 75.2%. Corrected lengths for 
other maps are provided in Additional file 2, Table S2. 

We observed little segregation distortion in F2 nauplii 
(Figure 2). Following Bonferroni correction for multiple 
tests, only one marker (TC078, on Chromosome 8) 
deviated significantly from Hardy- Weinberg equilibrium, 
with a deficiency of heterozygotes. We also observed 
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Table 1 Name, pre-amplification forward and reverse primers, extension primers and nucleotide calls for the 45 SNP 
loci used in this study 



Name 


Forward Primer 


Reverse Primer 


Extension sequence 


SD 


SC 


TC006 


GGGTATAGGGTTTGATCAAC 


TTGAGGACCTTTTACGCAAC 


AAGACCTGTGAAGAA 


G 


T 


TC008 


GGTTGATTGATTTTGGCTCC 


TGGAGTGCAAAAGTCCATCG 


GCCTGCATAGGCTCC 


C 


T 


TC011 


ATGAAAGACTCTTAAGTCG 


GATGTCTCACAAATCTCGGC 


ACTCTTAAGTCGACAGAC 


T 


A 


TC012 


GCAAACAGGACCGTGTTGAA 


ACAAGATGGAACGGGTAGAG 


CGTACACCTTGATCTTGCC 


C 


G 


TC016 


GTGAGI I I I G I CTCAAACGG 


ATTGCGAAGCCAAGTAAGCG 


TTTGTCTCAAACGGTTCACCCA 


G 


A 


TC017 


ATGCTTGGCAATTTGGAGTG 


TTTAGCTCCAGGACATTGGG 


GCAATTTGGAGTGTTCTGA 


C 


T 


TC033 


TCGCCAAGAGCAATAATCGG 


GATGGAAG I I I I GGTAGGTC 


TGAGAAAATCTCCTACTTACC 


C 


G 


TC040 


GAGGAAGTAGAGGCAATAGC 


ACCACGATGGAAGTTGCTTG 


CCAATAGCAACCATCA 


c 


T 


TC043 


GCATTGAGGATOTGGCTTC 


TGAGGAAAAATCACGTCCGC 


ATTGTGGCTTCACACCTA 


c 


T 


TC045 


ATTGGGATCCGACCCTGTTG 


GGTTTGAACACCTGAAATGG 


GGGATCTGGTTGTC 


T 


C 


TC046 


GGGAAAACGGGCTCAAGATG 


AAAGCTGCCCTTGTAAACGG 


ACAAATTGGCTTTGCTTGG 


C 


T 


TC051 


TCAAGCGAAACACCACCATC 


TCGTACAGCTGAATGGTGAC 


CTTCACCACCTACTC 


c 


T 


TC060 


GAGGAGCCA I I I I ICCACAC 


TCTTTTACAAGCACTCCTC 


GAAGTCATCACACACG 


G 


A 


TC073 


ACTTGGATGTGGTGTCCTTG 


CCAACACCAATGTCACCATC 


GTGGTGTCCTTGGCCTCA 


T 


A 


TC074 


TTGCTCAAGTCAACCATGCC 


GCAACCAGTTGATAGCTTCG 


TCCTGTArrrrGCCT 


C 


G 


TC077 


GTCTCAAACGCGAATGCAAG 


TTGCTTGGTTCGAGGACATC 


AGAAGGCACGGCCACG MUG 


T 


C 


TC078 


GGGTACTGTGGTGTCGATTC 


AAGAGGACAGCAAATGGGAC 


TGTGGATTCTCCTCGG 


C 


G 


TC084 


GCAAATTCTGGTCCACCATC 


ATGGAGGTTTATGGCTACCG 


TTTATCGGTACTGGTCTC 


T 


C 


TC085 


GAAGCCGACATCAAAATGAC 


AAGGAGAAATCTCGCATTGC 


CATCTGAACCAGAAATCGATC 


C 


A 


TC099 


TAGCAATCTTCGGTCACCAC 


CAGCCCAUATGATTGCTCG 


GGGTTTGCACGGGATA 


T 


C 


TC102 


GCGTCAAGGATGTTGAGAAC 


TGGTGCCCATACATGAGTTC 


GAACTTCATGGCCTTGAAGAT 


A 


T 


TC103 


GAGATCTCAGGAGAGTTGC 


AATGAGGAATCTTCTCTGGC 


GAGTTGCAACATCTC 


T 


C 


TC104 


GCTTACTTCATTCACGCACC 


GAGACTTTGAGCGTCTGAAC 


CCATTAmAGTCTATCCAAGA 


A 


G 


TC106 


GATTTCTTCTCCCACACCAG 


AATGAGAAGCAATGCTCCTG 


ACGGGGATCTTGAC 


C 


T 


TC107 


GCTTCCCATAAAGGTGAATC 


TGGAAAGAAACGGGAAGTGC 


AAAGGTCAATCAATGACT 


C 


A 


TC111 


AAGTGTCTCTGCGAGTGCTG 


TTCGAATGATGGGATCCGTG 


CAACCTTGGTCGGATGTTCCGGG 


c 


G 


TC112 


TGGGTGCCmGAAGAGAAC 


GAAACATCAAAGCAAGCAGG 


AGAGAACATCACCTAA 


c 


T 


TC118 


TCAAGGACACCATCAAGAGC 


TCCGGACAAGATGTTCTTGG 


AGAAGGACGCCTCCTA 


T 


C 


TC124 


ACACCAAACCATTGAGTGGC 


ATATGTAGAAGACCGAGAGC 


TACCCTGTTGATGCGA 


C 


T 


TC125 


TTGTTGGACrTGATGCCTCG 


GCTCGTATTGCTCGATCTTG 


ATGTGGACAAATTGAA 


G 


A 


TC128 


TGTCGTATCTCGTCCTCTrC 


AGAACACGAGAAGAAATTC 


CTCGTCCTGTTCTTGGAATC 


G 


A 


TC130 


ACGTTCGTTTGTTAAGGCGG 


GAACATTGCTGACAACGGAG 


GTTTGTTAAGGCGGTAACTG 


G 


A 


TC152 


GGCCCCAGACGACGACGA 


ATGCATTCGCGAAAGGCCTG 


CTCCCTCAAAGCCTC 


A 


G 


TC155 


TTGGATCAGATCCAGCGGT 


GCGAGGATTTGCAAGTACAC 


AGATCCAGCGGTGCCTTAAACC 


A 


C 


TC156 


AATTTGATCTAGGAGCAACG 


GCATTCAGTTCCAGTITGAG 


GGAGCAACGTTGCATTG 


C 


T 


TC157 


rAAGTTrCTGGrrAAGAAAG 


TrrATGAArATGTrrATGT 


GGGAACGATGATTTGGAA 


Q 


Q 


TC162 


GAAGATGCTTGAGTGGATGG 


TCGGCGGAGGTGTACACTIT 


AAGGCAATCTCCAGGAT 


T 


c 


TC167 


TGGGATCTTGATG I I I ICGC 


ACTGAGGAAGGATCTTGAGG 


AGGTTGGTCATTGT 


T 


c 


TC169 


GTAGCGCCTCAAGGTTATGG 


TCACGCTTGGCAACATGGAG 


AAGGGCCAATATGTGGC 


T 


c 


TC171 


ATGTCTCAAAACGGCCTTCC 


TCGAGGAAGTTTCGTTGCAC 


GCCTTCGGGATCAGGGGTA 


T 


c 


TC180 


AATGTCTCGCAAGAGCTGG 


TGGTGGTTACCTTTCGAAAT 


GAGCTGGACCCAGACAG 


T 


c 


TC184 


CTGCCGAATATCATGAAGGC 


GGGTCCTCAACGGAAATAAC 


AAGCCTCAGATTTGGGAGT 


G 


c 


TC188 


GGTCATTCTTCGAGACCTAC 


CGATAC I I I I I I GAGAGTAGG 


AGTGAACATGTCTGT 


G 


A 


TC189 


AGACGTTATGAAGCACGACC 


TTTGTTATCTTCCAAGGGC 


GCACGACCCAATCCGTGTC 


T 


A 


RPOL 


CAGGTTCA7TCAAAGCGGAC 


AATCCAAGATTTCAACCTC 


GGTTGGTAAGTCGG 


G 


T 



Abbreviations refer to the following T. californicus populations: SD: Point Loma, San Diego, California (32°45'N, 117°15'W; n = 28), SC: Santa Cruz, California 
(36°57'N, 122°03'W; n = 28). Additional file 1, Table S1 provides information on additional SNP loci and SNP calls for additional populations. 
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•10 c |v| I Linkage map from newly-hatched larvae (n=190) 

H Linkage map from adult males (n=205) 

Figure 1 Tigriopus californicus linkage maps generated from an F2 mapping cross between populations SD and SC The maps include 
45 SNPs and 8 microsatellites (1203, 56J2, 1555, 30, 228, 480, 558 and 197). Linkage groups with numerical names correspond to the numbered 
chromosomes identified in [31]. Maps generated from nauplii and adult males are both shown. 



some heterozygote deficiency at the linked marker 
TCI 56 and a deficiency of SC homozygotes at two mar- 
kers on Chromosome 1. The observed deficiency of het- 
erozygotes at microsatellite 1555 may be a technical 
artifact, due to a large peak size difference between the 
SC and SD alleles, in combination with generally lower 
peaks in nauplii compared to adults, causing heterozy- 
gotes to be under-called. Overall in nauplii we observed 
significantly more SD alleles than expected with a 1:1 
ratio assuming equal genetic contributions from both 
parental lines (j 2 = 21.3, p < 0.001). 

In contrast, we observed a high level of significant seg- 
regation distortion in adult males, distributed non-ran- 
domly across the genome and biased in different 
directions on different chromosomes (Figure 2). Follow- 
ing Bonferroni correction, eleven markers deviated sig- 
nificantly from Hardy- Weinberg equilibrium. All 
markers on Chromosome 10 exhibited a very strong 



deficiency of SD homozygotes. There was also a defi- 
ciency of SD homozygotes, associated with an excess of 
SC homozygotes, at TC011 on Chromosome 9. In con- 
trast, markers on Chromosomes 2 and 7 exhibited an 
excess of heterozygotes, associated with a deficiency of 
SC homozygotes. We also observed a deficiency of SC 
homozygotes and an excess of SD homozygotes at Chro- 
mosome B, although this was not significant at the Bon- 
ferroni corrected p value. In adults, we observed overall 
fewer SD alleles than expected with a 1:1 ratio (x = 
21.3, p < 0.001). 

We observed four deviations from expected two-locus 
multiplicative genotypic frequencies between physically 
unlinked loci at p < 0.05. Although it is unclear at what 
level to correct for multiple testing in this case [45], as 
comparisons involving different physically linked mar- 
kers are not independent, we note that none of these 
deviations were significant after a Bonferroni correction 
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□ SD/SC DSD/SD DSC/SC 

Figure 2 Proportions of heterozygotes, SD homozygotes and SC homozygotes observed at each marker locus for (i) newly hatched 
nauplii and (ii) adult males. The cytoplasmic background for all individuals was SD. Significance of deviations from expected 2:1:1 Mendelian 
segregation ratios is indicated by: *** p < 0.00094 (Bonferroni corrected p = 0.05); ** p < 0.01; * p < 0.05. Markers are arranged by linkage 
group, as indicated under chart. 



based on the number of possible comparisons between 
linkage groups (132). Three of these deviations involved 
marker TC167 on Chromosome C and markers on 
Chromosome 7 (TC 162, TC 184 and RPOL, Figure 3). 
Both individuals homozygous for SD and those homozy- 
gous for SC at Chromosome 7 markers exhibited a 
strong deficiency of SD homozygotes and a less marked 
excess of SC homozygotes at marker TC167. Those het- 
erozygous at Chromosome 7 markers showed an excess 
of SD homozygotes and a deficiency of SC homozygotes 
at TC167. The additional significant deviation was 
between marker TC012 on Chromosome D and TC017 
on Chromosome 2. Individuals homozygous for SC at 
TC012 exhibited a deficiency of SD homozygotes and 
an excess of heterozygotes at TC017; those homozygous 
for SD at TC012 showed an excess of SD homozygotes 
and a deficiency of heterozygotes at TC017 = 17.7, p 
= 0.02). 



In the non-recombinant backcross, several chromo- 
somes exhibited deviations from the expected 1:1 ratio 
of homozygotes (SCSC) to heterozygotes (SCSD), how- 
ever only one (Chromosome 10) was significant follow- 
ing Bonferroni correction, probably reflecting the small 
sample size (Figure 4). Chromosome 10 exhibited a dra- 
matic deficiency of heterozygotes; in contrast, heterozy- 
gote excess was observed at Chromosome 7 and 
Chromosome B. We found no evidence for epistatic 
interactions between chromosomes. 

Discussion 

In this study, we mapped 45 newly developed SNP mar- 
kers, and eight pre-existing microsatellite markers onto 
the genome of T. californicus. The map includes mar- 
kers on all 12 T. californicus chromosomes and has an 
estimated genomic coverage of 75%. This is the first 
linkage map to be generated for the Copepoda, one of 



Pritchard et at. BMC Genetics 201 1, 12:54 
http://www.biomedcentral.eom/1 471 -2 1 56/1 2/54 



Page 9 of 1 3 




-6 J 1 

■ TC1 84 _H_TC1 03 i □ RPOL 0 TC162 □ TC051 

Figure 3 Deviations from expected two-locus genotypic frequencies for marker TC1 67 and the five markers on Chromosome 7 The y 

axis indicates the contribution of each possible genotypic combination, as shown above the chart, to the total % 2 value. Negative value indicate 
that a genotypic combination was observed less frequently than expected by chance, while positive values indicate that a genotypic 
combination was observed more frequently. Total y} values were as follows: TC167 - TC184, j} = 17.1, p = 0.03; TC167 - TC103, y} = 14.5, p = 
0.07;TC167 - RPOL, % 2 = 16.2, p = 0.04; TC167 -TC162, % 2 = 18.3, p = 0.02; TC167 - TC051,x 2 = 11.9, p = 0.16. 



the most abundant animal taxa on the planet [46], and 
one of few currently available for crustaceans. All SNP 
markers are within coding loci, which will facilitate 
future investigations of chromosome synteny between T. 
californicus and other taxa. 

We used this map to investigate patterns of segrega- 
tion distortion across the genome of F2 hybrids between 
two T. californicus populations. We observed some 




0 I™i™i M i™i™i™i™i™i M i™i™i™ 

125789 10 11 ABCD 

Chromosome 



Figure 4 Proportions of heterozygotes (SDSC) and 
homozygotes (SCSC) observed at each chromosome in adult 
males from the non-recombinant backcross with cytoplasmic 
background SD. Significance of deviations from expected 1:1 
Mendelian segregation ratios are indicated by *** p < 0.004 
(Bonferroni corrected p = 0.05); ** p = 0.01. 



segregation distortion in newly hatched larvae, at fewer 
than 10% of markers at p < 0.05, which may be caused 
by pre-zygotic effects or pre-hatching mortality, and was 
not maintained until adulthood. Segregation distortion, 
however, was strikingly higher in adult males, affecting 
45% of marker loci at p < 0.05. This is comparable to 
the amount of segregation distortion that has been 
observed within interpopulation or inter-species crosses 
in other taxa (e.g. Nasonia spp, 29% of markers in adult 
males [45]; Mimulus guttatus 48% of markers in mature 
plants [47]; Arabidopsis lyrata, 50% of markers [48]; 
Daphnia magna, 33% of markers [44]; Daphnia pulex, 
21% of markers [49]; Lepomis spp. 36.8% of markers in 
fry [50]). As seen in other studies, distorted markers 
were not distributed randomly across the genome but 
clustered together on particular linkage groups. In T. 
californicus, this clustering may have been exacerbated 
by the lack of recombination in females, which increases 
the likelihood that markers physically linked to regions 
selected against in hybrids will be inherited together. 

Our observation of high segregation distortion in F2 
hybrid adults but not in F2 nauplii mirrors the pattern 
seen for several coding loci in different T. californicus 
interpopulation crosses (ME1, ME2 [29]; CYC, RISP, 
CYC1 [30,41]). This suggests that, in general, observed 
deviations from Hardy- Weinberg equilibrium in interpo- 
pulation T. californicus crosses are the result of selection 
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against genotypes between hatching and adulthood, 
rather than being due to meiotic drive or differential 
gametic fertilization success. Although segregation dis- 
tortion in inter-species and interpopulation hybrids is a 
common finding, our study is one of rather few to expli- 
citly demonstrate a role of postzygotic selection in gen- 
erating this phenomenon. Launey and Hedgecock [51] 
showed that segregation distortion in the oyster Crassos- 
trea gigas is caused by post-hatching mortality of indivi- 
duals homozygous for deleterious recessives. Martin and 
colleagues [52] identified hybrid genetic combinations 
affecting survivorship in Iris. Rogers and Bernatchez [53] 
found evidence for selection against particular genotypic 
combinations acting between fertilization and hatching 
in backcross hybrids between lake whitefish (Coregonus 
clupeaformis) ecotypes. Niehuis and colleagues [45] 
found evidence that cytonuclear co-adaptation caused 
genotypic-specific mortality between hatching and adult- 
hood in F2 interspecific hybrids in the wasp Nasonia. 

While we have data for only one of the two reciprocal 
crosses, our T. californicus results show very little evi- 
dence for cytonuclear coadaptation. All F2 individuals 
had an SD mitochondrial background: while we did 
observe an overall excess of SD alleles in nauplii, this 
had become an overall deficiency of SD alleles in adult- 
hood. In adults we found only two markers (on Chro- 
mosome B) where the homozygote matching the SD 
homozygote was favored, and eight markers where the 
pattern of segregation distortion was opposite to that 
which would be expected if the nuclear and mitochon- 
drial genomes within the two populations were co- 
adapted. For example, we observed dramatic segregation 
distortion in adult F2 males throughout Chromosome 
10, with the direction of segregation distortion indicat- 
ing strong selection against the nuclear genotype, SDSD, 
that matched the mitochondrial background. We also 
observed a large excess of SCSC homozygotes at Chro- 
mosome 10 in the non-recombinant backcross, suggest- 
ing that this genotype is more fit than the alternative, 
SDSC. Harrison and Edmands [31] correspondingly 
observed a deficiency of SDSD Chromosome 10 homo- 
zygotes in males (but not in females) in backcrosses 
between SD and another population, RP (Royal Palms). 
Taken together, these results suggest that part or all of 
Chromosome 10 derived from the SD population has a 
deleterious effect on viability in males of this interpopu- 
lation cross, that appears to act in an incompletely 
dominant manner. There are several reasons why such 
an apparently deleterious portion of the genome may be 
maintained in the SD population. First, this deleterious 
effect may be only expressed in a hybrid nuclear genetic 
background. Our study did not detect epistatic interac- 
tions involving Chromosome 10, although this may be 
due to limited power. Second, as we only examined 



males, it is possible that aspects of SD Chromosome 10 
may be advantageous in females. Alternatively, genes on 
SC Chromosome 10 may cause masculinization; we 
note, however, that while F2 offspring in this study did 
not deviate from a 1:1 sex ratio, the concurrently gener- 
ated nonrecombinant backcross was significantly female 
biased despite all individuals containing at least one 
copy of the SC chromosome (V.L. Pritchard and cowor- 
kers, unpublished data). Third, the deleterious aspect of 
Chromosome 10 may not be expressed in the natural 
environment. It may be masked in the wild SD popula- 
tion by the presence of a more dominant allele that was 
lost both from our isofemale lines and the SD parental 
line used in [31]. Even if this is not the case, previous 
studies with T. californicus have shown varying experi- 
mental conditions to alter the viability of different 
hybrid genotypes [26,27,54]. Additionally, the outcome 
of replicated experimental hybridizations may vary even 
under apparently identical conditions [41], suggesting 
that even apparently minor environmental changes can 
have a large influence on the fitness of different geno- 
types. Finally, T. californicus populations in the wild 
experience repeated population bottlenecks, which are 
expected to affect the outcome of selection. Hence even 
if the deleterious aspect of Chromosome 10 is expressed 
in the natural environment it may persist in the wild SD 
population due to drift. Indeed, there is evidence that 
many T. californicus populations carry such a genetic 
load [28]. For the SD population in particular, previous 
studies have indicated a selective disadvantage to SD 
homozygotes for the coding loci ME2 and CYC, even on 
the SD mitochondrial background [29,41]. In contrast, 
SD homozygotes were favored on a mismatching mito- 
chondrial background, for the coding loci ME1 and 
RISP [29,41]. We note, however, that these results vary 
by sex and study, and none of these four coding loci are 
on linkage groups exhibiting significant segregation dis- 
tortion in the current cross (ME1, Chromosome C; 
ME2, Chromosome A; CYC, Chromosome D; RISP, 
Chromosome 8, Rose and Edmands, unpublished data). 
Studies of taxa other than T. californicus have also 
shown that homozygotes mismatching the cytoplasmic 
background can be favored in hybrids. Fishman and col- 
leagues [55], for example, examining segregation distor- 
tion in Mimulus F2 hybrids, found a strong excess of M. 
guttatus homozygotes on a M. nasutus cytoplasmic 
background. Similarly Martin and colleagues [52], in a 
backcross study using Iris brevicaulis and /. fulva, found 
that, at three QTLs, presence of /. fulva homozygotes 
decreased long-term survivorship despite a matching 
cytoplasmic background. 

In contrast to the pattern observed for Chromosome 
10, we observe heterozygote excess, with no apparent 
selection against SDSD homozygotes, throughout most 
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of Chromosome 2 and Chromosome 7. We also observe 
an excess of heterozygotes for Chromosome 7 in the 
nonrecombinant backcross. This chromosome contains 
the locus coding for mtRPOL, which has been the focus 
of recent studies of cytonuclear coadaptation in T. cali- 
fornicus. In comparison to our results, Ellison and Bur- 
ton [22], looking at allelic frequencies in F4 hybrid 
adults, found evidence for selection against the SD 
mtRPOL genotype in crosses with both SC and another 
population, AB. They also observed that, unlike for 
other crosses, recombinant inbred lines with matching 
SD mtRPOL and SD mtDNA did not demonstrate the 
same OXPHOS transcriptional profile under conditions 
of hypo-osmotic stress as SD parentals, suggesting that 
an epistatically interacting nuclear locus is involved in 
mitochondrial transcription in SD. They suggested the 
transcription factor TFAM as a possible candidate. In 
this context, it is interesting that our results are sugges- 
tive of an epistatic interaction between marker TC167 
and Chromosome 7; recent work (Rose and Edmands, 
unpublished) has revealed marker TC167 to be closely 
linked to TFAM. Nevertheless, the genotypic association 
patterns between mtRPOL and TC167 are not what 
would be expected if there is simple co-adaptation 
between SD mtRPOL and SD TFAM; individuals homo- 
zygous for SD mtRPOL exhibit a deficiency of SD 
homozygotes, and a slight excess of SC homozygotes, at 
TC167. 

Higher divergence between parental lines is expected to 
result in increased frequencies of distorted loci [56], par- 
ticularly due to heterozygote deficits [49]. For example, 
parental divergence is cited as the reason why crosses 
between mildly divergent D. pulex populations result in 
21% of markers showing segregation distortion, largely 
due to homozygote deficits [49], while crosses between 
highly divergent Daphnia magna populations result in 
33% of markers showing transmission ratio distortion, 
largely due to heterozygote deficits [44]. This pattern is 
consistent with the prediction that overdominance 
between alleles of closely related taxa may yield to under- 
dominance between alleles in more distantly-related taxa 
[57]. Alternatively, divergence may increase the ratio of 
epistatic interactions involving heterozygous loci. In the 
current study hybridization between highly differentiated 
populations (over 20% mitochondrial divergence [15,16]) 
led to a high frequency of marker distortion in adults 
(45%), but no significant heterozygote deficits, indicating 
relatively slow accumulation of underdominance and/or 
epistasis involving heterozygotes. 

As has been seen in other Tigriopus studies [30], and 
in other taxa [58] it is clear that nuclear loci can inter- 
act in a complex way to influence fitness; unfortunately 
we lack the power to investigate such interactions in 
more depth in the current study. Additionally, we did 



not consider epigenetic effects, which have previously 
been suggested to alter gene transcription in interpo- 
pulation hybrids of T. californicus [59]. Overall, our 
results suggest many intriguing avenues for further 
investigation into the genetic basis of reduced fitness 
in interpopulation hybrids of T. californicus. These stu- 
dies will be greatly facilitated by the recent transcrip- 
tome assembly for both the SD and SC T. californicus 
populations [60], and by continuing advances in crus- 
tacean genomics [61]. 

Conclusion 

We developed 45 population-diagnostic SNP markers 
for Tigriopus californicus, with which we generated the 
first linkage map available for the Copepoda. We used 
this to examine segregation distortion in F2 interpopula- 
tion hybrids, which are known to have reduced fitness 
compared to parental populations. We found dramatic 
segregation distortion in adult males, but not in newly 
hatched larvae, indicating that this distortion arises as a 
result of selection against particular genotypic combina- 
tions between hatching and adulthood. Distorted mar- 
kers were not distributed randomly across the genome 
but clustered within particular linkage groups. In con- 
trast to other studies, we found little evidence for cyto- 
nuclear co-adaptation in this interpopulation cross. 
Instead, different linkage groups exhibited markedly dif- 
ferent patterns of distortion, that appear to have been 
influenced by nuclear-nuclear epistatic interactions and 
may also reflect genetic load carried within the parental 
lines. 

Additional material 



Additional file 1: Table SI: Name, accession numbers, primers, SNP 
calls and Blast matches for 51 nuclear and 3 mitochondrial loci 

Abbreviations refer to the following J. californicus populations: SD: Point 
Loma, San Diego, California (32°45'N, 117°15'W, n = 28); SC: Santa Cruz, 
California (36°57'N, 122°03'W, n = 28); PBJ: Punta Baja, Baja California 
(29°58'N, 1 15°48'W, n = 31); SUN: Sunset Beach, Washington (48°30'N, 
123°00'W, n = 20); FHL: Friday Harbor Laboratories, Washington (48°33'N, 
125°08'W, n = 15); LC: Leo Carillo Beach, California (34°03'N, 1 18°56 W, n 
= 20); PD: Point Dume, California (34°00'N, 118°48'W, n = 18); CAT: 
Catalina Island, California (33°27'N, 118 29'W, n = 19); SCI: Santa Cruz 
Island, California (34°03'N, 119°34'W, n = 17); RP: Royal Palms, Palos 
Verdes, California (33°42'N, 118°19'W, n = 13); AB: Abalone Cove, Palos 
Verdes, California (34°44'N, 118°19'W, n = 20); LB: Laguna Beach, 
California (33°33'N, 117°47'W, n = 19); OSP: Osprey Point, San Diego, 
California (32°44N, 117 °15'W, n = 16); PBN: Punta Banda, Baja California 
(31°43'N, 116°43'W, n = 19); PMO: Punta Morro, Baja California (31°52'N, 
1 16°40'W, n = 20); PA: Playa Altamira, Baja California (28°33'N, 114°05'W, 
n = 20). An 'n' indicates successful SNP calls in less than 80% of the 
individuals screened. Loci marked with * exhibit a deficiency of identified 
heterozygotes in known SD x SC Fl individuals (under-calling). 'Genbank 
# (primer design)' refers to the accessions from which sequencing 
primers were designed. 'BLAST match' shows results of a BlastX search of 
the NCBI protein database using these sequences; only matches with an 
E-value < 1 x F 5 are shown. 'Genbank # (sequence)' references the 
sequences from which the SNPs were identified. 
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Additional file 2: Table S2: Number of inter-locus comparisions (n), 
map distance in cM Kosambi (Distance), standard error of map 
distance (SE), and LOD scores for all markers on the twelve 
chromosomes (Chrom). Data are shown for nauplii and adults, both 
including and omitting microsatellite markers. Total map lengths and 
map lengths corrected following [43,44] are also provided. 



List of abbreviations 

ATP: adenosine triphosphate: cDNA: complementary deoxyribonucleic acid: 
CYC: cytochrome c; dNTP: deoxyribonucleotide triphosphate; LOD: logarithm 
of the odds; ME: malic enzyme; mRNA: messenger ribonucleic acid; PCR: 
polymerase chain reaction; QTL: quantitative trait locus; RISP: rieske iron- 
sulfur protein. 
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